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ABSTRACT 

The  yield  line  analysis  method  is  used  extensively  throughout  the  blast  design  and  analysis 
community  to  determine  the  strength  of  reinforced  concrete  structural  elements  to  resist  blast 
overpressures.  Technical  literature  has  long  been  available  that  describes  the  yield  line 
method  of  analysis  and  explains  how  to  derive  the  necessary  analysis  equations  for  a 
structural  element  such  as  a  wall  or  a  roof.  Most  presentations  of  the  yield  line  analysis 
method  deal  with  common  slab  configurations,  such  as  one-way  spans  of  varying  support 
conditions  or  two-way  spans  that  usually  consist  of  a  rectangular  slab  supported  along  two  or 
more  edges.  Very  little  discussion  is  available  on  how  to  analyze  slabs  with  openings,  and 
even  less  information  is  available  on  how  to  analyze  slabs  with  covered  openings.  Yet  slabs 
with  covered  openings,  such  as  blast-resistant  doors  or  windows,  are  very  common  structural 
elements  in  the  explosive  safety  design  community.  This  paper  presents  a  method  by  which 
the  yield  line  analysis  method  can  be  used  to  analyze  and  design  slabs  and  plates  with  open¬ 
ings. 

EFFECT  OF  OPENINGS  ON  STRENGTH  OF  SLABS 

In  blast  resistant  applications,  the  size  of  openings  in  slabs  can  become  very  significant  when 
determining  ultimate  resistance  to  pressure.  Openings  tend  to  attract  yield  lines,  but  they 
don't  automatically  weaken  a  slab.  Openings  result  in  less  surface  area  to  collect  load,  and 
under  some  conditions,  the  ultimate  strength  of  a  slab  can  actually  increase.  When  a  slab  has 
openings  that  are  covered  there  is  no  reduction  in  surface  area  collecting  load  and  the  designer 
has  to  add  into  the  analysis  the  effects  of  the  additional  blast  load  collected  by  the  cover. 

These  cover  loads  are  passed  into  the  slab  via  line  loads  around  the  supported  edges  of  the 
cover.  These  line  loads,  if  the  cover  is  large  enough  or  strong  enough,  can  dramatically 
reduce  the  strength  or  ultimate  resistance  of  the  slab.  Yield  line  analysis  methods  that  include 
the  effects  of  line  loads  can  more  accurately  determine  the  strength  of  a  uniform  slab  and 
often  eliminate  the  need  for  pilasters  and  headers  around  the  openings. 
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GENERAL  PRINCIPLES  OF  YIELD  LINE  ANALYSIS 


Yield  line  analysis  is  founded  upon  the  principle  of  conservation  of  energy:  the  work 
performed  by  an  external  force  moving  through  a  distance  is  equal  to  the  internal  work 
performed  by  rotations  about  plastic  hinges  that  resist  the  external  force.  The  yield  line 
analysis  method  provides  an  upper  limit  estimate  of  the  maximum  ultimate  resistance,  r^,  of  a 
slab  for  an  assumed  mode  of  failure.  The  slab  is  assumed  to  fail  by  deflecting  until  plastic 
hinges  form  along  the  supports  and  in  the  interior  of  the  slab.  These  plastic  hinges  or  yield 
lines  subdivide  the  slab  into  planar  sectors  that  rotate  about  these  hinges  until  maximum 
deflection  or  failure  occurs.  Any  of  several  geometric  combinations  of  plastic  hinges  can 
describe  a  valid  mode  for  failure  (failure  mechanism)  of  a  given  slab,  and  each  must  be 
checked  to  determine  which  describes  the  lowest  value  of  r^.  The  failure  mechanism  that 
defines  the  lowest  value  of  r^  will  require  the  least  amount  of  energy  to  fail.  The  assumption 
of  an  incorrect  failure  mechanism  for  a  slab  will  result  in  either  overpredicting  the  ultimate 
resistance  of  the  element  or  in  an  undefined  solution.  This  underscores  the  importance  of 
checking  all  of  the  credible  failure  mechanisms  for  a  slab  to  ensure  the  most  reliable  r^  has 
been  identified. 

There  are  two  methods  available  in  yield  line  analysis  that  can  be  used  to  calculate  r^:  the 
virtual  work  method  and  the  equilibrium  method.  Each  method  has  its  advantages  and 
disadvantages  as  a  design/analysis  tool. 

VIRTUAL  WORK  METHOD 

The  virtual  work  method  of  yield  line  analysis  makes  use  of  the  principle  "if  a  rigid  body  that 
is  in  equilibrium  under  a  system  of  forces,  is  given  a  virtual  displacement,  the  sum  of  the 
virtual  work  done  by  the  forces  is  zero"  (Ref  1).  A  failure  mechanism  consisting  of  plastic 
hinges  (yield  lines)  and  planar  sectors  is  assumed  for  a  given  slab.  When  an  external  pressure 
(blast  load)  is  applied  to  the  slab  and  it  begins  to  deflect,  yield  lines  begin  to  form.  As  the 
yield  lines  form,  they  become  the  axis  of  rotation  of  the  planar  sectors  formed  out  of  the  slab. 
Each  sector  is  assumed  to  act  as  a  rigid  body.  As  each  sector  rotates  and  deflects  under  the 
pressure,  external  and  internal  work  is  performed. 

External  Work 


External  work  is  represented  by  the  external  forces  acting  on  the  slab  moving  through  the 
distance  the  slab  deflects.  The  amount  of  external  work  performed  is  best  calculated  by 
computing  the  external  forces  acting  on  each  rigid  sector  and  multiplying  them  by  the 
displacement  of  that  sector.  This  is  done  by  multiplying  the  pressure  acting  on  the  slab  by  the 
area  of  the  sector,  then  multiplying  that  product  by  the  distance  the  centroid  of  the  area  of  the 
sector  has  displaced. 

When  a  slab  contains  a  covered  opening,  the  cover  over  the  opening  resists  the  external 
overpressure  and  transfers  that  pressure  to  the  slab  via  a  line  load  along  the  cover  supports. 
The  external  work  performed  by  the  line  loads  along  a  opening  is  computed  by  multiplying 


the  load  times  the  length  of  the  support  to  define  a  force.  The  force  is  then  multiplied  by  the 
distance  the  line  load  moved.  That  distance  can  accurately  be  described  as  the  deflection  of 
the  slab  at  the  centroid  of  the  line  load. 


The  external  work  performed  over  the  entire  slab,  becomes  the  sum  of  the  external  work 
performed  on  each  sector  by  the  overpressure  plus  the  sum  of  the  work  performed  by  any  line 
loads  around  openings  in  the  slab: 


EQUATION 


where: 


Af  =  Area  of  Sector  i 
6^  =  Deflection  at  the  centroid  of  the  sector 
V  =  Line  load  per  unit  length 
Lj  -  Length  of  the  line  load  j 
6,  =  Deflection  at  the  centroid  of  line  load  j 


To  determine  how  much  overpressure  load  the  slab  can  safely  resist,  the  pressure  acting  on 
the  slab  is  defined  as  the  ultimate  resistance,  r^,  of  the  slab. 

Internal  Work 


As  the  sectors  form  and  rotate,  the  ultimate  moment  capacity  of  the  slab  resists  the  bending 
moments,  shear,  and  torsional  forces  along  the  yield  lines.  This  is  defined  as  internal  work 
and  can  be  computed  by  multiplying  the  ultimate  moment  capacity  of  the  slab  along  each 
yield  line  (the  ultimate  moment  capacity  per  unit  length  times  the  length  of  the  yield  line)  by 
the  angle  of  rotation  along  the  yield  line.  The  total  internal  work  for  a  slab,  is  defined  as 
the  sum  of  the  moment  capacity  along  all  yield  lines: 


EQUATION 


(2a) 


where:  M  =  Ultimate  moment  capacity  along  yield  line 
/■  =  Length  of  yield  line  i 
0^  —  Rotation  about  yield  line  i 

The  angle  of  rotation,  0,  can  be  approximated  for  small  rotations  by  tan  0.  In  most 
cases  where  steel  reinforcement  is  parallel  and  perpendicular  to  the  supports,  the  internal  work 
equation  becomes 


EQUATION 


(2b) 


where: 


Mtc 

^vy 


=  Ultimate  moment  capacity  along  yield  line  i  in  x  direction 
=  Length  of  yield  line  i  projected  in  the  y  direction 
=  Rotation  about  yield  line  i  projected  about  the  y  axis 
=  Ultimate  moment  capacity  along  yield  line  i  in  y  direction 
=  Length  of  yield  line  i  projected  in  the  x  direction 
=  Rotation  about  yield  line  i  projected  about  the  x  axis 


The  total  external  work  is  set  equal  to  the  total  internal  work  for  a  slab.  The  resulting 
equation  can  then  be  solved  for  r^.  An  assumed  failure  mechanism  is  considered  to  be  invalid 
if  the  solution  to  the  virtual  work  equation  yields  a  negative  value. 

Advantages  and  Disadvantages  of  Virtual  Work 

The  advantage  of  the  virtual  work  method  is  that  it  treats  the  entire  slab  as  a  unit  which 
simplifies  somewhat  the  derivation  of  the  equations.  Imbalances  of  shear  and  torsion  forces 
along  the  yield  lines  cancel  themselves  out  when  forces  are  summed  across  the  entire  slab. 

The  disadvantage  of  the  virtual  work  method  is  that  the  exact  locations  of  the  yield  lines  are 

not  known  and  they  must  be  solved  for.  The  equation  derived  to  describe  r^  for  slabs  with 
openings  is  large,  complex,  and  usually  nonlinear.  The  exact  solution  to  the  equation  requires 


partial  differentiation  with  respect  to  each  of  the  unknown  locations  of  yield  lines.  The  result¬ 
ing  equations  are  large,  cumbersome  and  error  prone.  A  trial  and  error  solution  where  the 
absolute  minimum  r^  is  computed  based  on  assumed  locations  of  the  yield  lines  is  an 
alternative  solution  method. 

EQUILIBRIUM  METHOD 

The  equilibrium  method  is  similar  to  the  virtual  work  method,  except  that  the  external  and 
internal  work  are  computed  and  equated  on  a  sector  by  sector  basis.  Imbalances  of  shear  and 
torsional  forces  along  the  boundaries  between  sectors  exist  and  must  be  accounted  for  in  the 
equilibrium  method.  The  adjustments  for  the  imbalance  of  forces  between  sectors  are  known 
as  nodal  forces.  Nodal  forces  in  rectangular  slabs  with  uniform  reinforcement  in  each 
direction  are  significant  only  when  yield  lines  intersect  a  free  edge  or  an  opening  such  as  a 
door  or  window  (Figure  1). 


Figure  1.  Nodal  Forces. 


Figure  1.  Nodal  Forces. 


The  expression  for  the  nodal  force,  Q,  is 


EQUATION 

Q  =  ±Mj^cot0 

where:  =  Moment  capacity  of  steel 

perpendicular  to  free  edge 
0  =  Acute  angle  formed  by 

intersection  of  yield  line  with 
free  edge 


Note  that  when  the  yield  line  intersects  a  free  edge  at  a  90-degree  angle,  the  nodal  force  is 
zero.  The  nodal  force  is  negative  in  the  acute  angle  and  positive  in  the  obtuse  angle. 

An  expression  for  r^  is  derived  for  each  sector  under  this  method.  A  solution  is  reached  when 
the  locations  of  the  yield  lines  are  adjusted  such  that  the  value  for  r^  is  the  same  for  each  of 
the  sectors.  The  solution  can  be  solved  directly  by  equating  the  expressions  for  r^  from  each 
sector  and  solving  them  as  simultaneous  equations.  This  effort  is  complex,  as  each 
relationship  for  r^  is,  as  with  the  virtual  work  method,  usually  a  complex,  nonlinear,  algebraic 
expression. 

An  assumed  failure  mechanism  is  considered  invalid  when  a  simultaneous  solution  of  the 
equilibrium  equations  yield  a  negative  value  for  r^  or  a  yield  line  location.  The  failure 
mechanism  is  also  invalid  when  the  equations  will  not  converge  on  a  common  value  for  r^ 
when  using  a  trial  and  error  solution. 

COMPARISON  OF  VIRTUAL  WORK  AND  EQUILIBRIUM  METHODS 

Mathematical  equations  derived  to  solve  for  r^  of  any  particular  failure  mechanism 
theoretically  yield  basically  the  same  result  whether  derived  using  the  equilibrium  method  or 
the  virtual  work  method.  To  assure  accuracy  and  minimize  potential  errors,  the  equations  for 
each  failure  mechanism  described  in  this  paper  were  derived  and  checked  using  both 
methods.  Sample  problems  were  run  to  verify  that  solutions  from  either  method  gave 
essentially  the  same  result.  In  every  failure  case  identified  the  results  for  r^  agreed  within  1% 
or  less.  The  locations  of  the  yield  lines  at  times,  however,  did  vary  significantly.  In  cases 
where  differences  did  occur,  the  location  of  yield  lines  predicted  by  the  equilibrium  method 
were  assumed  to  be  the  most  correct  because  the  work  energy  is  balanced  equally  among  all 
sectors. 

The  virtual  work  equation  is  somewhat  simpler  to  derive,  but  is  more  complex  to  solve 


directly.  The  equilibrium  equations  can  be  solved  simultaneously,  but  involve  more  complex 
algebraic  equations. 

A  trial  and  error  solution,  though  it  can  be  used  for  both  methods,  works  best  with  the 
equilibrium  method.  As  yield  line  locations  are  assumed,  the  computed  value  of  r^  for  each 
sector  identifies  trends  that  indicate  which  direction  the  assumed  locations  must  be  adjusted  to 
achieve  convergence.  The  logical  identification  of  these  "trends"  work  nicely  with  computer 
programming  algorithms.  The  trial  and  error  procedure  is  repeated  until  the  designer  is 
satisfied  that  the  optimum  locations  of  yield  lines  have  been  identified  that  yield  a  minimum 
value  for  r^. 

Because  of  the  complex  nature  of  the  equations,  yield  line  analysis  solutions  of  slabs  with 
covered  openings  are  only  practical  when  used  in  the  form  of  a  computer  program. 

REDUCTION  OF  MOMENT  CAPACITY  IN  CORNERS 

Yield  lines  first  begin  to  form  in  a  structural  element  along  the  supports  and  on  the  interior  of 
the  slab.  As  deflection  increases  and  a  collapse  mechanism  begins  to  form,  the  yield  lines 
extend  and  grow.  Yield  lines  that  form  on  the  slab  interior  (positive  yield  lines)  eventually 
extend  to  intersect  with  those  that  form  along  the  supports  (negative  yield  lines).  At  this  point 
the  collapse  mechanism  has  fully  formed.  Experimental  test  results  over  the  years  clearly 
show  that  as  yield  lines  extend  into  the  corners  of  a  slab,  the  localized  stiffness  prevents  full 
rotations  from  developing,  and  the  moment  capacity  of  the  slab  in  its  corners  is  not  fully 
developed.  This  causes  the  yield  line  analysis  to  overestimate  r^  in  these  cases.  Some  design 
criteria  recommend  that  reduction  factors  be  applied  to  the  computed  r^  to  more  closely 
predict  the  actual  value  of  the  ultimate  resistance  (Ref  1,  2).  The  magnitude  of  the  reduction 
factor  varies  according  to  the  type  of  reinforcement  and  the  geometry  of  the  slab.  The  DOD 
design  criterion  (Ref  3)  recommends  that  the  moment  capacity  be  reduced  by  1/3  in  corner 
regions  to  account  for  this  reduction. 

When  a  1/3  reduction  of  moment  capacity  is  accounted  for  in  the  derivation  of  yield  line 
equations  for  slabs  with  openings,  already  complex  equations  become  significantly  more 
complex.  For  this  reason,  the  yield  line  equations  presented  here  were  first  derived  without 
considering  any  reduction  of  slab  moment  capacity  in  the  corners.  Once  agreement  on  r^  was 
reached  between  the  two  methods,  the  equations  were  expanded  to  account  for  the  decreased 
capacity  to  absorb  energy  in  the  corners.  Several  important  observations  were  noted  when  the 
1/3  reduction  of  moment  capacity  was  added  to  the  equations.  (1)  The  solution  for  r^  was 
consistently  predicted  to  be  10-15%  lower  than  that  predicted  without  the  reduction  in 
moment  capacity.  (2)  The  solution  for  r^  using  the  virtual  work  method  would  not  converge 
exactly  with  the  equilibrium  method.  The  virtual  work  equation  consistently  yielded  a  1-2% 
lower  value  for  r^.  (3)  The  virtual  work  method  predicted  significantly  different  locations  for 
the  yield  lines.  These  phenomena  were  observed  even  in  rectangular  slabs  without  openings. 
The  locations  of  the  yield  lines  predicted  by  the  equilibrium  method  were  the  same  with  or 
without  a  reduction  of  moment  capacity  in  corners.  The  reason  for  these  differences  have  not 
been  determined,  but  they  are  most  likely  related  to  an  inadvertent  introduction  of  an 


imbalance  of  equilibrium  forces  between  internal  and  external  work. 


To  minimize  differences  between  solutions  predicted  by  the  two  methods,  and  to  simplify 
complex  equations  as  much  as  possible,  a  singular  reduction  in  r^  to  account  for  corner  effects 
is  recommended  as  discussed  in  References  1  and  2.  Table  1  shows  recommended  factors  for 
rectangular  slabs  with  right  angle  corners. 


Table  1.  Recommended  Reductions  in  Ultimate  Resistance  of  Slab,  r^. 


Type  of 

Slab 

Reduction  per  Corner  (%)  | 

No  Restraining  Moments  on  Edges 

Restrained  Edges, 
Top  Reinforced 
in  Comers 

Top  Reinforced 
in  Corners^ 

No  Top 
Reinforcement 

All  Sides  Supported 

0 

-1 

-2 

One  Long  Free  Edge 

-1 

-3 

-3 

“Top  steel  area  per  unit  width  equal  to  bottom  steel  area  per  unit  width. 


Table  1.  Recommended  Reductions  in  Ultimate  Resistance  of  Slab,  r„. 


YIELD  LINE  CONFIGURATIONS 

The  most  likely  failure  mechanisms  for  slabs  with  openings  along  with  the  equations  that 
have  been  derived  for  those  mechanisms  are  given  in  Appendix  A.  For  any  slab  with  an 
opening,  these  failure  patterns  at  a  minimum  should  be  investigated  as  possible  collapse 
mechanisms.  Simple  geometrical  considerations  will  eliminate  most  failure  mechanisms  for 
any  given  slab.  Other  possible  mechanisms  may  exist  and  must  be  determined  by  following 
yield  line  theory  and  using  good  engineering  judgment. 

The  configurations  shown  in  Appendix  A  are  designed  to  take  advantage  of  symmetry  for 
economy  of  space.  If  a  slab  does  not  exactly  match  the  geometry  shown  in  the  figures, 
chances  are  it  will  match  one  of  the  cases  if  the  orientation  is  rotated  90  degrees  and/or  re¬ 
versed.  For  example,  a  slab  with  a  door  in  the  lower  left  corner  can  be  reversed  to  match  the 
opening  in  the  lower  right  corner  shown  in  Cases  1  through  4.  If  an  opening  in  a  slab  is  closer 
to  a  side  wall  than  the  floor,  the  failure  pattern  for  the  opening  shown  in  Case  1 8  may  apply 
by  rotating  the  slab  90  degrees  and  adjusting  the  input  to  the  equations  to  match  the  new  rotat¬ 
ed  position.  Similarly,  slabs  with  openings  should  be  rotated  as  many  times  as  necessary  to 
match  up  with  as  many  of  the  failure  mechanisms  as  possible  to  ensure  that  all  possible  cases 
have  been  checked.  Software  can  easily  be  designed  to  perform  these  checks  automatically. 


The  yield  line  equations  shown  in  Appendix  A  assume  that  the  slab  has  a  covered  opening.  If 
a  slab  is  to  be  analyzed  that  has  no  cover  over  the  opening,  the  line  load  forces  from  the  cover 
are  simply  input  as  zero.  This  will  cause  all  of  the  line  load  terms  to  drop  out. 

Modes  of  Failure 


A  slab  with  a  covered  opening  can  fail  in  one  of  two  modes.  The  first  and  most  common 
mode  of  failure  is  when  the  slab  fails  because  of  overpressure.  The  second  mode  of  failure 
occurs  when  the  slab  fails  primarily  due  to  line  loads  from  a  cover  over  an  opening.  This 
mode  of  failure  is  not  as  common,  but  is  important  for  the  designer  to  be  aware  of. 

A  rigorous  yield  line  analysis  of  a  slab  with  a  covered  opening  will  then  consist  of  two 
phases.  The  first  phase  will  assume  a  failure  mechanism  based  on  overloading  due  to 
overpressure.  The  ultimate  resistance  for  the  slab  is  then  computed,  accounting  for  the  effects 
of  blast  overpressure  and  line  loads  acting  on  the  slab.  The  second  phase  assumes  a  yield  line 
failure  mechanism  based  on  an  overload  from  line  loads  acting  on  the  wall.  The  ultimate 
resistance  for  the  slab  is  then  recalculated  accounting  for  both  the  blast  overpressure  and  the 
line  loads  on  the  slab  based  on  that  failure  mechanism.  The  lowest  value  of  r^  controls  the 
design. 

This  paper  documents  only  the  first  phase  which  represents  the  most  common  situations 
encountered  in  design  and  analysis  of  slabs  with  openings.  To  determine  how  much  line 
loads  influence  the  ultimate  resistance,  analyze  the  slab  with  and  without  line  loads. 

If  the  r^  with  line  loads  is  less  than  50%  of  r^  without  line  loads,  a  phase  2  analysis  should 
also  be  considered. 

DERIVATION  OF  YIELD  LINE  EQUATIONS 

To  perform  a  yield  line  analysis,  a  failure  mechanism  is  first  assumed  for  a  slab.  This  paper 
shows  the  derivations  of  the  equations  used  to  describe  the  failure  mechanism  shown  in 
Figure  2  (See  Case  2  in  Appendix  A).  This  is  a  reinforced  concrete  slab  of  height,  H,  and 
length,  L,  that  is  fixed  along  all  four  edges,  with  a  door  of  height,  a,  and  width,  b,  located  in 
the  lower  right  corner.  The  door,  while  resisting  a  blast  pressure  will  apply  line  loads, 
across  the  top  support,  and  across  the  side  support.  Yield  lines  are  identified  by  dashed 
lines.  The  variables  w,  y,  and  z  define  the  locations  where  yield  lines  intersect  free  edges  of 
the  doorway.  The  wall  slab  is  divided  by  the  yield  lines  into  four  sectors  identified  1-4. 


Figure  2.  Yield  Line  Failure  Mechanism  for  Case  2. 


Figure  2.  Yield  Line  Failure 
Mechanism  for  Case  2. 


The  yield  line  equations  derived  here  equate  the  external  work  with  the  internal  work  as  the 
structure  deforms  under  a  blast  load.  The  sectors  defined  are  assumed  to  be  rigid  bodies  with 
planar  rotation  about  its  supports.  Yield  lines  that  form  along  the  supports  are  called  negative 
yield  lines,  and  those  forming  in  the  span  between  sectors  are  called  positive  yield  lines.  The 
greatest  displacement  obtained  by  any  one  sector  of  the  slab  is  defined  as  a  unit  displacement. 
All  other  displacements  throughout  the  slab  are  normalized  and  expressed  in  terms  of  the  unit 
displacement.  If  an  opening  is  located  on  a  portion  of  the  slab  where  a  maximum 
displacement  would  otherwise  occur,  the  yield  lines  are  extended  into  the  opening  until  they 
intersect.  The  hypothetical  intersection  becomes  the  assumed  location  of  the  maximum, 
hence,  the  unit  displacement. 


Figure  3.  Extension  of  Yield  Lines  to  Intersection. 


Figure  3 .  Extension  of  Yield 
Lines  to  Intersection. 


Figure  4.  Horizontal  Extension  of  Yield  Lines  from  Point 


Figure  4.  Horizontal  Extension  of  Yield 
Lines  From  Point  e. 


For  example,  in  Figure  3,  Line  a-b  separates  Sectors  1  and  2.  Line  c-d  separates  Sectors  2 
and  3.  These  two  lines  are  extended  to  intersect  each  other  at  Point  e.  The  maximum  or  unit 
displacement  is  assumed  to  occur  at  Point  e.  Sectors  3  and  4  are  separated  by  Line  f-g.  Line  f- 
g  must  now  be  extended  to  intersect  the  positive  yield  Line  c-e. 

There  are  two  possibilities  of  intersection.  The  first  possibility  is  Line  f-g  will  extend  directly 
to  intersect  somewhere  along  Line  c-e.  To  maintain  planar  rotation  of  Sector  3,  Line  f-g  must 
be  constrained  to  only  intersect  Line  c-e  between  Points  d  and  e.  The  other  possibility  is  that 
the  extension  of  Line  f-g  does  not  intersect  Line  c-e  at  all,  but  passes  somewhere  to  the  right 
of  it.  In  this  case  a  horizontal  line  must  be  extended  out  to  the  right  from  Point  e.  Line  f-g  is 
then  extended  to  intersect  the  horizontal  line.  This  defines  Point  h  as  shown  in  Figure  4.  The 
displacement  along  Line  e-h  is  defined  as  the  unit  displacement,  and  the  rotations  of  all 
sectors  as  planar  sections  remain  consistent. 

Preliminary  geometrical  relationships  that  are  needed  in  the  derivation  of  the  yield  line 
equations  are  defined  as  illustrated  in  Figure  5  and  as  given  below: 
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Figure  5. 

Geometric  Relationships  Used  in  Derivation  of  Equations  for  Case  2. 


(a)  ^2  +  W2  >  L.  (b)  ^2  ^2 

Figure  5.  Geometric  Relationships  Used  in  Derivation  of  Equations  for  Case  2. 


For  Sector  1,  the  equation  for  external  work  was  derived  using  Equation  1  and  is  shown  in 
Table  2(a).  The  internal  work  equation  was  derived  using  Equation  2a  and  is  shown  in  Table 
2(b).  Setting  and  solving  for  yields  the  relationship  shown  in  Table  2(c).  The 

units  used  in  these  equation  derivations  are  feet  for  distances  and  lengths,  inch-pounds  per 
inch  for  moments,  psi  for  ultimate  resistance,  and  pounds  per  foot  for  line  loads.  A 
conversion  factor  of  144  is  used  in  all  external  work  terms  to  keep  units  consistent. 

Eor  Sector  2,  the  equation  for  external  work  is  shown  in  Table  3(a),  while  the  internal  work 
equation  is  shown  in  Table  3(b).  Setting  solving  for  r2  yielded  the  equation 

shown  in  Table  3(c). 

Two  equations  are  possible  for  Sectors  3  and  4,  depending  on  the  location  of  the  yield  lines. 

If  X2  +  W2  >  L  as  shown  in  Eigure  5(a),  the  equations  are  derived  as  follows. 

Eor  Sector  3,  the  equation  for  external  work  is  shown  in  Table  4(a),  while  the  internal  work  is 
shown  in  Table  4(b).  Setting  solving  for  is  shown  in  Table  4(c). 

Eor  Sector  4,  the  external  work  is  shown  in  Table  5(a),  while  the  internal  work  is  shown  in 
Table  5(b).  Setting  and  solving  for  yields  the  equation  shown  in  Table  5(c). 


If  X2  +  W2  <  L,  the  equations  for  Sectors  3  and  4  are  derived  from  the  geometry  shown  in 
Figure  5(b).  The  equations  are  derived  as  follows. 

For  Sector  3,  the  equation  for  external  work  is  the  same  as  shown  in  Table  4(a).  However, 
the  internal  work  is  different  as  shown  in  Table  4(e).  Setting  and  solving  for  is 

shown  in  Table  4(f). 

For  Sector  4,  the  external  work  equation  is  shown  in  Table  5(d),  while  the  internal  work 
equation  is  shown  in  Table  5(e).  Setting  and  solving  for  yields  Table  5(f). 

The  relationships  derived  for  r^,  r2,  r3,  and  are  all  assumed  to  be  equivalent  (r^  =  r2  =  r3  =  r4). 
These  equations  may  now  be  solved  simultaneously  or  by  trial  and  error  methods. 


_ Table  2.  Yield  Line  Equations  for  Sector  1 

(a)  External  Work: 


(b)  Internal  Work: 
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Table  2.  Yield  Line  Equations  for  Sector  1 


Table  3.  Yield  Line  Equations  for  Sector  2 


Table  3.  Yield  Line  Equations  for  Sector  2 


Table  4.  Yield  Line  Equations  for  Sector  3. 
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Table  4.  Yield  Line  Equations  for  Sector  3 


Table  5.  Yield  Line  Equations  for  Sector  4. 
If  ^2  +  W2  >  L,  then 


(a)  External  Work: 


=  144rA±UH-a)w 
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(b)  Internal  Work: 


(c)  Ultimate  Resistance: 
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If  ^2  +  ^>2  ^2  L,  then 


(d)  External  Work: 


(e)  Internal  Work: 


Table  5.  Yield  Line  Equations  for  Sector  4. 

Equations  were  derived  for  Cases  1,3-21  using  similar  procedures  and  are  listed  in  Appendix 
A. 


STIFFNESS  OF  SLABS  WITH  OPENINGS 


The  stiffness  of  slabs  with  openings  can  have  different  levels  of  importance  when  computing 
the  maximum  deflection  caused  by  a  blast  load.  Stiffness  is  most  important  when  expected 
maximum  slab  deflections  will  be  in  the  elastic  and  elasto-plastic  ranges.  Stiffness  becomes 
less  important  when  maximum  deflections  extend  more  into  the  plastic  deformation  range 
where  the  ultimate  resistance,  r^,  increasingly  influences  the  response.  When  maximum 
plastic  deflections  are  large,  the  value  for  stiffness  can  be  approximated  with  little  loss  of 
accuracy  by  computing  the  stiffness  of  a  slab  or  plate  as  though  there  is  no  opening. 

When  greater  accuracy  is  required  in  computing  the  stiffness  of  a  slab  with  an  opening,  the 
following  approach  is  recommended.  Once  a  yield  line  analysis  has  been  performed,  a  failure 
pattern  identified,  and  the  location  of  the  yield  lines  determined,  two  equivalent  slabs  are 
created.  The  first  equivalent  slab  (Slab  1)  will  have  the  same  height,  width,  and  support 
conditions  as  the  original  slab.  Slab  1  will  have  a  uniform  moment  capacity  in  the  horizontal 
and  vertical  directions  and  for  positive  and  negative  bending.  The  uniform  moment  capacity 
is  determined  by  averaging  M^ni,  The  stiffness  of  this  slab  is 

computed  using  standard  procedures  and  neglecting  the  effect  of  openings. 

Case  1  -  Case  4 


For  Cases  1-4  shown  Appendix  A,  compute  the  ratios  a/H  and  b/L.  If  a/H  is  greater  than 
b/L,  then  compute  the  stiffness  of  a  second  equivalent  slab  (Slab  2)  fixed  across  the  top, 
bottom  and  left  side,  with  a  free  edge  along  the  right  side  as  shown  in  Figure  6.  Use  a  height, 
H,  and  a  length  L-b.  Use  the  uniform  moment  capacity  for  the  equivalent  slab  computed 
above.  Compute  the  difference  between  stiffness  of  the  two  slabs  by  subtracting  the  stiffness 
of  Slab  2  from  the  stiffness  of  Slab  1.  Multiply  this  difference  by  the  ratio  a/H  and  subtract 
the  result  from  the  stiffness  of  Slab  1. 

If  a/H  is  less  than  b/L,  then  compute  the  stiffness  of  Slab  2  assuming  the  slab  is  fixed  across 
the  top,  right  and  left  sides,  with  a  free  edge  along  the  bottom  as  shown  in  Figure  7.  Use  a 
height  of  H-a  and  a  length  L.  Use  the  uniform  moment  capacity  for  Slab  1  computed  above. 
Compute  the  difference  between  the  stiffness  of  the  two  slabs  by  subtracting  the  stiffness  of 
Slab  2  from  Slab  1.  Multiply  this  difference  by  the  ratio  b/L  and  subtract  the  result  from  the 
stiffness  of  Slab  1. 

Case  5  -  Case  12 


For  Cases  5-12  shown  in  Appendix  A,  compute  the  ratios  a/H,  c/L,  and  (L-b-c)/L.  If  a/H  is 
greater  than  c/L  and  (L-b-c)/L,  then  compute  the  stiffness  of  equivalent  Slab  2  assuming  it  is 
fixed  across  the  top,  bottom  and  left  side,  with  a  free  edge  along  the  right  side  as  shown  in 
Figure  6.  Use  a  height  H  and  a  length  c  or  L-b-c,  whichever  is  larger.  Use  the  uniform 
moment  capacity  for  the  equivalent  slab  computed  above.  Compute  the  difference  between 
stiffness  of  the  two  slabs  by  subtracting  the  Slab  2  stiffness  from  Slab  1. 

Multiply  this  difference  by  the  ratio  a/H  and  subtract  the  result  from  Slab  1 . 


Figure  6. 

Equivalent  Slab  2  -  Top,  Bottom,  and  Left  Sides  Fixed,  Right  Edge  Free. 


Figure  6.  Equivalent  Slab  2  -  Top, 
Bottom,  and  Left  Sides  Fixed, 
Right  Edge  Free. 


If  a/H  is  less  than  c/L  or  (L-b-c)/L,  then  compute  the  stiffness  of  equivalent  Slab  2  fixed 
across  the  top,  right  and  left  sides,  with  a  free  edge  along  the  bottom  as  shown  in  Figure  7. 
Use  a  height  of  H-a  and  a  length  L.  Use  the  uniform  moment  capacity  for  Slab  1  computed 
above.  Compute  the  difference  between  stiffness  of  the  two  slabs  by  subtracting  the  Slab  2 
stiffness  from  the  Slab  1  stiffness.  Multiply  this  difference  by  the  greater  of  the  two  ratios  c/L 
or  (L-b-c)/L  and  subtract  the  result  from  the  stiffness  of  Slab  1. 


Figure  7.  Equivalent  Slab  2  -  Top,  and  Sides  Edges  Fixed,  Bottom  Edge 

Free. 


Figure  7.  Equivalent  Slab  2  -  Top 
and  Side  Edges  Fixed,  Bottom 
Edge  Free. 


Case  13  -  Case  21 


For  cases  13-21  shown  in  Appendix  A,  compute  the  ratios  d/H,  (H-a-d)/H,  c/L,  and  (L-b- 
c)/L.  Ratio  1  is  the  larger  of  the  two  ratios  d/H  and  (H-a-d)/H.  Ratio  2  is  the  larger  of  the  two 
ratios  c/L  and  (L-b-c)/L. 

If  Ratio  1  is  less  than  Ratio  2,  then  compute  the  stiffness  of  equivalent  Slab  2  assuming  it  is 
fixed  across  the  top,  bottom  and  left  side,  with  a  free  edge  along  the  right  side  as  shown  in 
Figure  6.  Use  a  height  H  and  a  length  c  or  L-b-c,  whichever  is  larger.  Use  the  uniform 
moment  capacity  for  Slab  1  computed  above.  Compute  the  difference  between  stiffness  of 
the  two  slabs  by  subtracting  the  Slab  2  stiffness  from  the  Slab  1  stiffness.  Multiply  this 
difference  by  the  ratio  a/H.  Subtract  the  result  from  the  stiffness  of  Slab  1. 

If  Ratio  1  is  greater  than  Ratio  2,  then  compute  the  stiffness  of  equivalent  Slab  2  fixed  across 
the  top,  right  and  left  sides,  with  a  free  edge  along  the  bottom  as  shown  in  Figure  7.  Use  a 
height  of  d  or  H-a-d,  whichever  is  greater,  and  a  length  L.  Use  the  uniform  moment  capacity 
for  Slab  1  computed  above.  Compute  the  difference  between  the  two  stiffnesses  by 
subtracting  the  Slab  2  stiffness  from  the  Slab  1  stiffness.  Multiply  this  difference  by  the  ratio 
b/L.  Subtract  the  result  from  the  stiffness  of  Slab  1. 


SUMMARY 


New  design  aids  for  yield  line  analysis  of  slabs  with  covered  openings  have  been  presented. 
Because  of  the  complex  nature  of  yield  line  equations,  the  analysis  is  best  handled  via 
computer  programs.  FORTRAN  subroutines  for  each  of  the  failure  mechanisms  presented  in 
this  paper  are  currently  under  development  at  the  Naval  Facilities  Engineering  Service  Center. 
This  work  is  being  performed  as  part  of  the  Defense  Nuclear  Agency  DAHS-CWE  project. 
The  yield  line  equations  shown  here  were  derived  using  the  equilibrium  method.  Eor  the 
software  that  is  under  development,  the  equilibrium  method  is  being  used  to  predict  the 
locations  of  the  yield  lines.  The  virtual  work  method  is  also  being  used  to  compute  the 
ultimate  resistance  of  the  slab  as  it  tends  to  be  slightly  more  conservative. 

SAMPLE  PROBLEM 

Compute  the  ultimate  resistance  and  stiffness  of  a  rectangular  slab  with  a  door  as  shown  in 
Eigure  8.  The  wall  is  36  inches  thick  and  contains  #8  steel  reinforcing  bars  at  12  inches  on 
center  in  the  vertical  direction  and  #9  steel  reinforcing  bars  at  12  inches  on  center  in  the 
horizontal  direction.  The  door  is  constructed  of  a  2-inch-thick  A-36  steel  plate  supported 
across  the  top  and  two  sides  and  free  spanning  across  the  bottom.  The  wall  will  be  allowed 
up  to  8  degrees  of  support  rotation  due  to  blast  loading. 

Figure  8.  Wall  with  Covered  Door  Opening,  Sample  Problem. 


Figure  8.  Wall  With  Covered 
Door  Opening,  Sample  Problem. 


SOLUTION: 


STEP  1 .  Compute  the  moment  capacity  of  the  walls 


The  moment  capacity  is  calculated  using  standard  procedures  outlined  in  Reference  3. 


=  187,393  in.-lbs/in. 
Mh,2  =  187,393  in.-lbs/in. 
Mhp  =  187,393  in.-lbs/in. 


=  238,141  in.-lbs/in 
M,,2  =  238,141  in.-lbs/in 
M,p  =238,141  in.-lbs/in, 


STEP  2.  Compute  the  ultimate  resistance  of  the  door  and  corresponding  line  loads  on  the 
wall. 


From  Reference  3,  r^  of  the  door  is  computed  to  be  72.0  psi.  The  intersection  of  the  yield 
lines  occur  at  x  =  3.75  feet,  y  =  5.1 1  feet.  From  these  tributary  areas. 


EQUATIONS 


^10' +  5.11'^ 


(3.750(72p5/*12*12) 


10^ 


=  29,375  lbs//t 


7.5' (10' -5.1 1')  (72^5/ *12  *12)  *  25,345 


Ibs/ft 


STEP  3.  Compute  the  ultimate  resistance,  r^.  of  the  slab. 

Solve  for  r^  of  the  slab  for  each  of  Cases  1-4  in  Appendix  A.  Use  the  smallest  value  of  r^ 
obtained  from  the  four  cases.  The  results  of  a  trial  and  error  solution  for  Cases  1-4  are 
summarized  below: 


TABLE  6 


Case  Number 

Ultimate  Resistance, 

1 

189.1  psi 

2 

did  not  converge 

3 

did  not  converge 

4 

did  not  converge 

STEP  4.  Compute  the  average  moment  capacity  of  equivalent  Slab  1. 

EQUATION 

^  3(18L3g3).^141)  .  212,767 

STEP  5.  Compute  the  stiffness  of  Slab  1. 

Prom  procedures  in  Reference  3:  =  2053.9  psi-in. 

STEP  6.  Compute  the  ratios  a/H.  b/L:  Determine  configuration  of  Slab  2. 

a/H  =  10  ft/15  ft  =  0.66667  ft  b/L  =  7.5  ft/25  ft  =  0.30000  ft 

a/H  is  greater  than  b/L.  Therefore,  Slab  2  will  be  fixed  along  top,  bottom,  and  left  side,  with  a 
height  of  15  feet  and  a  width  of  17.5  feet.  =  212,767  as  before. 

STEP  7.  Compute  the  stiffness  of  Slab  2. 

Prom  procedures  in  Reference  3:  =  1223.5  psi-in. 

STEP  8.  Compute  the  difference  between  the  two  stiffnesses. 


K,  (Slab  1)  -  K,  (Slab  2)  =  2053.9  -  1223.5  =  830.4  psi-in. 


STEP  9.  Multiply  the  difference  by  the  ratio  a/H. 


830.4  (a/H)  =  830.4  (0.66667)  =  553.6  psi-in. 
STEP  10.  Subtract  the  result  from  the  stiffness  of  Slab  1. 


(Slab  with  opening)  =  189.1  psi 
Kg  (Slab  with  opening)  =  2053.9  -  553.6  =  1500.3  psi-in. 
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TABLE  OF  ABBREVIATIONS 


Geometric  Terms 


H 

L 

a 

b 

c 

d 


Height  of  the  slab,  ft. 

Width  of  the  slab,  ft. 

Height  of  opening,  ft. 

Width  of  opening,  ft. 

Distance  from  opening  to  right  side  wall,  ft. 
Height  of  opening  above  floor,  ft. 


j:’  o:’  S  S:  S  S:  S  S 


Structural  Properties 

=  Negative  moment  capacity  along  horizontal  support  No.  1,  in.-lb/in. 

=  Negative  moment  capacity  along  horizontal  support  No.  2,  in.-lb/in. 

=  Positive  horizontal  moment  capacity  within  the  slab,  in.-lb/in. 

=  Negative  moment  capacity  along  vertical  support  No.  1,  in.-lb/in. 

=  Negative  moment  capacity  along  vertical  support  No.  2,  in.-lb/in. 

=  Positive  vertical  moment  capacity  within  the  slab,  in.-lb/in. 

=  Ultimate  resistance  of  a  slab,  psi. 

=  Ultimate  resistance  of  Sector  1  of  a  slab,  psi. 

=  Ultimate  resistance  of  Sector  2  of  a  slab,  psi. 

=  Ultimate  resistance  of  Sector  3  of  a  slab,  psi. 

=  Ultimate  resistance  of  Sector  4  of  a  slab,  psi. 

=  Ultimate  resistance  of  Sector  5  of  a  slab,  psi. 

=  Horizontal  line  loads  along  either  side  of  opening,  Ib/ft. 

=  Vertical  line  loads  along  top  and/or  bottom  of  opening,  Ib/ft. 


Yield  Line  Terms 

u  =  Location  of  yield  line,  ft. 

V  =  Location  of  yield  line,  ft. 

w  =  Location  of  yield  line,  ft. 

X  =  Location  of  yield  line,  ft. 

y  =  Location  of  yield  line,  ft. 

=  Location  of  yield  line,  ft. 


z 


APPENDIX  A 


YIELD  LINE  FAILURE  MECHANISMS  AND  EQUATIONS 

CASES  1  THROUGH  21 


24  x(H -y)  +  72  (L -X  -  w)(H -y)  +  24  w(H -y) 


24Hw 


VIEWGRAPHS  1  THROUGH  27 
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Internal 


Problem  Parameters 


Wall  Thickness  =  36  in. 

Steel  Reinforcement  =  #8  bars  @  12  in.  o.c. 

#9  bars  @  12  in.  o.c. 
Allowable  support  rotation  =  8  deg 


Door  =  A-36  steel,  2  in.  thick 
Supported  along  top  and  two  sides 


Step  1 :  Compute  the  Moment  Capacity  of  Walls 

Use  standard  computer  program 
such  as  BARCS 

.  =  187,393  in-lbs/in 

.  =  187,393  In-lbs/ln 

.  =  187,393  in-lbs/in 

.  =  238,141  in-lbs/in 

.  =  238,141  in-lbs/in 


238,141  in-lbs/in 


Step  2:  Compute  Ultimate  Resistance  of  Door 
and  Corresponding  Line  Loads 


From  TM  5-1300 
and  BARCS: 

r^  =  72.0  psi 

=  29,375  Ibs/ft 

=  25,345  Ibs/ft 


Step  3:  Compute  Ultimate  Resistance  of  Slab 


Case  1 


Case  2 


Approximate  Location  of  Yield  Lines 


Summary  of  Ultimate  Resistance, 


Case 


r 


u 


1 

2 

3 


189.1  psi 
did  not  converge 
did  not  converge 
did  not  converge 


4 


Step  4:  Compute  Average  Moment  Capacity  and 
Stiffness  of  Equivalent  Slab  1 

M..,  =  [M,,  + 

+  /  6 

=  [3(187,393)  +  3(238,141)]  /  6 

=  213,767  in.-lb/in. 

From  TM  5-1300;  Computer  Program  BARCS: 
=  2053.9  psi-in. 


qe 


Step  6:  Compute  Stiffness  of  Slab  With 
Covered  Opening 

Difference  between  K  ,  and  K 

el  ez 

-K^2  =  2053.9  -  1223.5 

=  830.4  psi-in. 

Ke  =  K  1  -  (K,,  -  (a/H) 

=  2053.9  -  830.4  (0.6667) 

=  1500.3  psi-in. 


